Model of the bremsstrahlung emission accompanying interactions between protons 
and nuclei from low up to intermediate energies: role of magnetic emission 
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A new model of the bremsstrahlung emission which accompanies proton decay and collisions of 
protons off nuclei in the energy region from the lowest up to intermediate, has been developed. This 
model includes spin formalism, potential approach for description of interaction between protons and 
nuclei, and operator of emission includes component of the magnetic emission (defined on the basis 
of Pauli equation) . In the problem of the bremsstrahlung during the proton decay in the first time 
a role of the magnetic emission is studied using such a model. For the studied Tm nucleus it has 
been studied the following: (1) How much does the magnetic emission change the full bremsstrahlung 
spectrum? (2) At which angle is the magnetic emission the most intensive relatively electric one? 
(3) Is there some space region where the magnetic emission increases strongly relatively electric 
one? (4) How intensive is the magnetic emission in the tunneling region? (5) Which values has the 
probability at its maximum and at zero energy limit of the emitted photons? It is demonstrated 
that the model is able to describe enough well experimental data of the bremsstrahlung emission 
which accompanies collisions of protons off the ^C, ^''Cu and ^°^Ag nuclei at the incident energy 
Tiab = 72 MeV (at the photon energy up to 60 MeV), the ^Be, ^^C and ^°*Pb nuclei at the incident 
energy Tiab = 140 MeV (at the photon energy up to 120 MeV). 

PACS numbers: 41.60.-m, 03.65.Xp, 23.50.+Z, 23.20.Js 
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I. INTRODUCTION 

According to theory of collisions of protons off nuclei, interactions between two nucleons play important role, which 
become leading at increasing of energy. By such a way, interaction between two nucleons (i.e. nucleon-nucleon, 
or two-nucleons interaction) is putted into the basis of relativistic models of collisions, with further application of 
formalism of Feynman's diagrams. But, on the other side, consideration of nucleus as medium allows to include space 
distribution of all nucleons into the model. Such a way takes into account non-locality of quantum mechanics, one of 
its fundamental aspects. Comparing these two different considerations, a question arises: what is more fundamental, 
interaction between different point-like nucleons of the studied nuclear system or quantum effects of non-locality in 
it? 

How important non-local effect in study of many-nucleons interactions? How much are they small? Results of 
[l| give some answer on this question: it was shown that fully quantum consideration of the boundary and initial 
conditions in theproblem of proton decay has essential influence on the calculated half-life (for example, half-lives 
calculated in [3t1Q[ can be changed up to 200 times after taking such conditions into account, while assumed error 
is only some percents in that models). This estimation indicates that non-local effects are not so small and their 
inclusion into calculations is able to essentially change results sometimes. 

Another aspect is collective motion. Models with nucleon-nucleon interaction should be the most accurate, if the 
collective effects caused by interactions between nucleons of the complete nuclear system were very small. However, 
we know that this is not so at low energies. One can assume that many-nucleons interactions disappear at increasing 
of energy of interacting nucleons. If to analyze bremsstrahlung emission, which accompanies collisions of protons 
off nuclei, then there are indications that two-nucleons interactions give the largest intensity of emission. But, we 
find that many-nucleons effects should arise at increasing of energy of the emitted photons^. We find confirmation 
about essential influence of many-nucleons interactions on the process of emission and all importance of its study in 
literature (for example, see J56| : in particular, two-nucleons approaches do not give positive explanation of nature of 
hard photons). 

Properties of the bremsstrahlung accompanying scattering of protons off nuclei have been studied enough well 



^ For example, in the problem of a-decay at increasing of energy of the emitted photon for obtaining stable value of the emission probability, 
it needs to continuously increase external boundary of space region of integration. In the task of fission this problem is essentially more 
difficuh (see [ij'l. 



(for example, see review [Tj], also [Tj] for emission in collisions between heavy ions). As a rule, as emitter of 
photons in nuclear system, both the nucleus as medium, and different nucleons in it were considered. The process of 
emission is studied as result of deacceleration of motion of nucleons in the averaged field of nucleus or in consequence 
of nucleon-nucleon collisions. At the same time, it was pointed out (for example, see [56') that properties of the 
nuclear bremsstrahlung emission accompanying nucleon-nucleus and nucleus-nucleus collisions (especially, in region 
of intermediate energies up to 150 MeV / nucleon) have been studied worst of all. This causes our interest in use of 
the optical model potentials 11] and folding potentials [12] for investigations of the bremsstrahlung emission, which 
accompanies interactions protons with nuclei. It could be interesting to obtain the model, which allows to describe 
the spectra in energy region from minimal up to intermediate. Possibility to take quantum non-local properties into 
account in description of such interactions reinforces our interest in such potential approach. 

However, in inves tig ations of the bremsstrahlung emission, which accompanies a-decay of nuclei [15l-l40l|. spontaneous 
fission of nuclei |4ll45l| , ternary fission of nuclei [52| , and also collisions of nucleons off nuclei |53l - l5q | , ions and nuclei off 



nuclei at non-relativistic energies |14{ , the emission caused by the magnetic moment of the fragment moving relatively 
the nucleus has not been taken into account. Such a way could be explained, if at such energies of the emitted photons 
the magnetic emission is enough small and it can be neglected in calculations (for example, see [50]). Microscopic 
models can provide a powerful formalism for study of many-nucleons interactions, where wave functions were obtained 
from a single-configuration resonanting group calculations. But, in particular, we see that magnetic emission was not 
included into such models, which were applied for description of the bremsstrahlung emission during scattering of 
protons on a-particles |57| 1 , a-particles on a-particles and light nuclei |58l . |59[ . 

The magnetic emission is connected with magnetic momentum and spin of the fragment, interacting with nucleus. 
Attempt to take such aspects into account leads to matrix form of equations of interactions (where two-component 
Pauli equation is the simplest) and many-component wave function of nuclear system (for example, see [60| , p. 32- 
35, 48-60). However, the magnetic component of emission and spin formalism are included in relativistic models of 
collisions of nucleons between themselves and with nuclei at intermediate energies (based on Dirac equation). Here, I 
should like to note two directions of intensive investigations: Refs. [6l|,[62l and Refs. (63l - l74| . However, main emphasis 
in these papers was made on construction of correct relativistic description of interaction between two nucleons in this 
task, where formalism was developed in momentum representation mainly. So, it could be interesting to obtain the 
model, combining spin formalism of interacting fragments of nuclear system (with inclusion of magnetic momentum) 
and potential approach for description of interaction between themselves. 

The problem of the bremsstrahlung during collisions of protons off nuclei and proton decay can be convenient in this 
investigation. In [T^ the problem of the bremsstrahlung during proton-decay was studied (see also [Ig|)- However, 
here the magnetic emission caused by the magnetic moment of proton was not taken into account (but spin-orbital 
component of potential was included and its influence on the spectrum was estimated). In order to clarify its role, a 
model with such aspect is needed. Main aim of this paper is construction of such a model. 

What interesting and new could this model give? We shall put such questions. How much does the magnetic 
emission change the full bremsstrahlung spectrum? At which angle is the magnetic emission the most intensive 
relatively electric one? Is there some space region where the magnetic emission increases strongly relatively electric 
one? How intensive is the magnetic emission in the tunneling region? Which values has the probability at its maximum 
and at zero energy limit of the emitted photons? We answer on such questions in this paper. 

II. MODEL 

A. Operator of emission of the bremsstrahlung photon 

Let us consider generalization of Pauli equation for A+l nucleons of the proton- nucleus system in laboratory frame 
(obtained starting from eq. (1.3.6) in [601, p. 33) 

i'h— = H^, iJ=^|^— (p, -^A,) +z,eA,,o-^^cr.rotA, l+l/(ri...r^+i), (1) 

Z— 1 ^ ^ 

where we use for any nucleon with number i (Hke to eq. (1.3.4) in f60'| for one-particle problem) 

X=7^^{v^- —A.) i>. (2) 

2miC V c / 

Here, $ = (x, V') is bispinor wave function of the proton-nucleus system, rrii and zi are mass and charge of nucleon 
with number z, Aj is component of potential of the electromagnetic field formed by this nucleon (describing possible 



bremsstrahlung emission of photon caused by this nucleon), a are Puah matrixes, A is mass number of nucleus, 
V{y\ . . -r^+i) is potential of (nuclear and Coulomb) interactions between all nucleons^. We pass to the center-of- 
masses frame, where we have distance r between center-of masses of proton and nucleus (for example, see Appendix A 
in Ref. [52], also [5g|). Then, one can represent this hamiltonian as H — Hq + W^ where W combines all items 
of electromagnetic field, which we define as operator of emission of the bremsstrahlung photon, and Hq is rest of 
hamiltonian without the emission of photons. Neglecting by relative motion of nucleons of nucleus in calculation of 
W, we find: 
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where Z^^l and m are effective charge and reduced mass of the proton-nucleus system, p is operator of momentum 
corresponding to r. Neglecting items at e^A^/c^ and Aq, the operator of emission in Coulomb gauge can be rewritten 
as 



W = -Zcff Ap - Zcff CT 

mc Zmc 



e / h 

rot A = — Zcff — Ap + — (7 ■ rot A 

mc \ 2 



Substituting the following form of the potential of electromagnetic field: 




we obtain 
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Here, e^") are unit vectors of polarization of the photon emitted (e^"''* = e*^")), kph is wave vector of the photon 
and Wph = fcphC — kph c. Vectors e'"' are perpendicular to kph in Coulomb calibration. We have two independent 
polarizations e*^^' and e"' for the photon with impulse kph {a = 1,2). One can develop formalism simpler in the 
system of units where h — 1 and c — 1, but we shall write constants h and c explicitly. Also we have properties: 
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Matrix element of emission 



Let us consider the matrix element in form: 



F: 



fi 



W k, 



rfir)WMr)dr, 



(8) 



where ipi{r) — \kij and V'/(r) = 1^/) are stationary wave functions of the proton- nucleus system in the initial i-state 
(i.e. state before emission of photon) and final /-state (i.e. state after emission of photon) which do not contain 
number of photons emitted. Substituting the operator of emission in form ^ into eq. ([S]) , we obtain: 



Ff, = {kf\W\h) = Z, 



e 2Trhc^ r \ 

'.S \ / i Pel + Pmag.l + Pmag.2 f , 

mc y Wph '^ J 



(9) 



^ According to [60| (see p. 32), the equation JTJ is working if energy Si of any nucleon with number i is close to its mass mi, i.e. 
\ei — rrii \ <^ mi (c = 1). From here one can obtain high energy limit for proton incident energy Sp ■C 2mp ~ 1. 86 GeV. By other words, 
inside energy region up to Sp the equation ((ij includes all relativistic properties, which Dirac equation gives us (with application of 
eq. l|2}). In particular, this limit is essentially higher intermediate energies for proton- nucleus collisions studied in this paper. 
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This definition for Ffi is in compHance with our previous formalism in [30|, |31|, I36l - l40l l49l l52l |75| . In particular, for 
square of matric element of emission we have (see eqs. (l)-(2) in |75j): 



|a/,|2 = 2TTT\Ff,\'^-6{wf-w,+Wpi,). 



(11) 



C. Wave function of nuclear system and summation over spinor states 

We shall define the wave function of the proton in field of the nucleus. We shall construct it in form of bilinear 
combination of eigenfunctions of orbital and spinor subsystems (as eq. (1.4.2) in [601, P- 42). However, we shall assume 
that it is not possible to fix experimentally states for selected M (eigenvalue of momentum operator J^). So, we shall 
be interesting in superposition over all states with different M and define the wave function so: 



^,i{v,s) = R{r) Y. Y. ci:::ijr,^^Yu^,)v,{si 



(12) 



where R{r) is radial scalar function (not dependent on m at the same I), n^ = r/r is unit vector directed along r, 
yim(nr) are spherical functions (we use definition (28,7)-(28,8), p. 119 in 77]), Cf^^,^ are Clebsh-Gordon coefficients, 
s is variable of spin, M — m + fi and I ~ j ± 1/2. For convenience of calculations we shall use spacial wave function as 

^i„^{r)^Ri{r)Yi,n{nr). (13) 

Spinor function 11^1(5) has two components w^i(s) and v^i^{s), which are eigenfunctions of spin operator Sz having 
eigenvalues cti and 172 (see (TTj], p. 247). So, we have: 



Action of operator of spin on the wave function is given by (see eq. (55,4) in [731, P- 248) 



(14) 



(15) 



and we have non-zero matrix elements: 
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(16) 



From eqs. ([TE)) (at s = 1/2, a ~ ±1/2) we calculate: 



(17) 



and find summations: 

Si,Sf=±l/2 Si,Sf=±l/2 Si,Sf=±l/2 



Considering vectorial form of spin operator, these formulas can be rewritten as 



s/ = ±l/2 



=±1/2 
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(19) 



.,-1/2 > +ez, 



Si,Sf=±l/2 



where orthogonal unit vectors Gx, Gy, g^ are used. 

So, using the found eqs. (J18p - (I19I) . we perform summation in eqs. (jlO[) over all spinor states: 
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where (fcj . . . fc^) is one-component matrix element 
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R} (r) YymA^r)* I R^ (r) Yi^mA^r) dv, 



(21) 



where integration should be performed over space coordinates only. 

We orient frame vectors Gx, Gy and e^ so, that g^ be directed along to kph. Then, vectors Gx and Gy can be directed 
along g'^' and g'^^ correspondingly. In Coulomb gauge we obtain: 

Gx = g(i), Gy = g(2), |ex| = |Gy| = |gz| = 1, |g(3)| = 0. (22) 

Now we perform summation in eqs. (|10p over polarization vectors and obtain: 
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D. Matric elements integrated over space coordinates 

We shall calculate the following matrix elements: 
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^*(i.)e-*^p.>'-^(p,(r)dr. (24) 
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1. Expansion of the vector potential A by multipoles 



Let us expand the vectorial potential A of electromagnetic field by multipole. According to [78| (see (2.106), p. 58), 
in the spherical symmetric approximation we have: 



C^e 
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where (see [78|, (2.73) in p. 49, (2.80) in p. 51) 
A;p„^(r,M) == jf^Jfcphr) T/p,jp,^^^(nr) 
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' ^ph . /, N rp . N (26) 
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Here, A; ,j;i(r, M) and A; j,;i(r, i?) are magnetic and electric multipoles, ji ,^ (fcphr) is spherical Bessel function of order 
Zph, T;^!^;' ,^(n].) are vector spherical harmonics. Eq. (P5| is solution of the wave equation of electromagnetic field in 

form of plane wave, which is presented as summation of the electrical and magnetic multipoles (for example, see p. 83- 
92 in 60*). Therefore, separate multipolar terms in eq. (j25l) are solutions of this wave equation for chosen numbers jph 
andlph (jph is quantum number characterizing eigenvalue of the full momentum operator, while Zph = Jph— 1, iph, jph+1 
is connected with orbital momentum operator, but it defines eigenvalues of photon parity and, so, it is quantum number 
also). 

We orient the frame so that axis z be directed along the vector kph (see [7g, (2.105) in p. 57). According to [7^ 
(see p. 45), the functions T; j^c .^(ni.) have the following form (^o = 0): 



Tjph;ph:™(ni-) = X! ('ph:l>Jph|™-M,Ai,"i) rip„,„_^(n,.) C/x, 



(27) 



where (/, 1, j rn — /i, /i, m) are Clebsh- Gordon coefficients, Yim{d, f) are spherical functions defined, according to [77 
(see p. 119, (28,7)-(28,8)). From eq. ([25]) one can obtain such a formula (at e"^^) = 0): 
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2. Spherically symmetric decay 
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Using dMl), for dMl) we find 
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Now we shall calculate components in eqs. ((23|) . For the first and third items we obtain: 
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Now we shall analyze the second item in eqs. p3| and find 
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Taking properties (jA7l) into account, we calculate eq. (p4|) further and obtain: 
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So, we have found all components in ([1 
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3. Calculations of the components pf'^^^, pf^^^ and pf^^^, pf^^^ 
For calculation of these components we shall use gradient formula (see [78|, (2.56) in p. 46): 
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and obtain: 
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^. Differential (angular) matrix elements of emission 
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We shall be interesting in the angular emission of the bremsstrahlung photons. By such a reason let us introduce 
e following differential matrix elements, dpf^ and dpf , dependent on the angle 9: 
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^ I {h + l){lph + l) dlE{k,lf,lph,h + l,lph- 1,A*) (jn 1 1 n i id i i ^\\ 

+ \/(2;, + i)(2V + i) ^^^ede [Mh.hJp^-i) -h-J2{h.if,ip^-i)]- 

~ V(2;, + i)(2;ph + i) S^^^d^ |Jift,^/,/ph + i) -k-Uh.if,ip^. + i)], 

(43) 
and 

"PjphM _ / ?ph + 1 dl {li,lf,lpi„lpi,~ l,^) ~ _n_ / /ph c^^(^i,^/,^ph,^ph + 1,m) t., , , -.X 

sin^d^ " Y2/ph + l sin^de •Jl'^,i/,(ph ij ]^ 2lpu + 1 sin0d0 • J l«„ (/, iph + ij. 

(44) 
One can see that integration of such a defined functions over the 6 angle inside region from to tt gives the fuh matrix 
elements pf^j^„ and pf ^ „ defined by eq. (|38l). and matrix elements P*^^ and pf j^„ defined by eq. (l40l) . 



E. Angular probability of emission of photon with impulse kph and polarization e'"'' 

I define the probability of transition of the system for time unit from the initial z-state into the final /-states, 
bein g in the given interval dvf, with emission of photon with possible impulses inside the given interval dvph, so (see 
Ref. [Zl, (42,5) § 42, p. 189; Ref. [zi], § 44, p. 191): 

I 12 

dW = ^^ ■d:^^2n |F/,|' d{wf - w^ + u-ph) • diy, 

d^fcph wl^^dwphdnph E,~Ef 

dv = dVf ■ dl/ph, dl^ph = ^^ = ^^^^p , W, -Wf^ = Wf^, 

where dt'ph and dvf are intervals defined for photon and particle in the final /-state, dJlph = d cos 0ph = 
sin 0ph rfflph c^'Pph , fcph = Wph/c. However, we have to take into account that in multipolar expansion (I25p for the 
vectorial potential of the electromagnetic field we oriented the frame so that axis z be directed along the vector kph . 
So, we have to do not use dilp^ in eq. (j45p . Ffi is integral over space with summation by quantum numbers of the 
system in the final /-state. Such procedure is averaging by these characteristics and, so, Ffi is independent on them. 
Interval dvf has only new characteristics and quantum numbers, by which integration and summation in Ffi was not 
performed. Integrating eq. (|45p over dwfi and substituting eq. ^, we find: 



dW ^ 



ZlffC^ hwph 



p{h,kf] 



dwr,y,. (46) 



Wph 



This is the probability of the photon emission with impulse kph (and with averaging by polarization e'"' ) where the 
integration over angles of the particle motion after the photon emission has already fulfilled. 

I define the following probability of emission of photon with momentum kph when after such emission the particle 
moves (or tunnels) along direction n^: differential probability concerning angle is such a function, definite integral 
of which over the angle 9 with limits from to tt equals to the total probability of the photon emission ( [46| ). Let us 
consider function: 

jfWiOfl_ _ Z!,^-!^f^^,^^,^^dpJ.,kf,Of)^^ ^^^^ 



dwph dcosOf 2-Kc^ nn? \ dcoa6f 

This probability is inversely proportional to normalized volume V. With a purpose to have the probability independent 
on V, I divide eq. (|47p on flux j of outgoing a-particles, which is inversely proportional to this volume V also. Using 
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quantum field theory approach (where v{p) — |p|/po at c = 1, see [801, § 21.4, p. 174): 

h c^ ki 



j = niv{pi), V, = |v.,| = 






E^ 



(48) 



where n^ is average number of particles in time unit before photon emission (we have rii — 1 for the normalized wave 
function in the initial z-state), v{pi) is module of velocity of outgoing particle in the frame system where colliding 
center is not moved, I obtain the differential absolute probability (let us name dW as the relative probability): 



dP{ipf,9f) d^Wi^f,0f) 



E, 



dw 



'ph 



"e//' 

dwph d cos 0f h c^ ki 27r c^ 



ZL.e^ WphE, 



I ki 



P{ki,kf) 



dp*{h,kf,Q,f) 



d cos 6 



/ 



(49) 



Note that alternative theoretical way for calculations of the angular bremsstrahlung probabilities in a-decays was 
developed in [35] based on different definition of the angular probability, different connection of the matrix element 
with the angle 9 between fragment and photon emitted, application of some approximations. 
Using formula pO|) . we rewrite eq. p9)) as 



(50) 



dP{ipf,0f) _ dPcl{ipf,9f) dPmag,liVf,9f) dPmag,2(',5/,^/) , d PintcifcrcnceiV> f , 9 f) 



dWph 



dwph 



dwph 



dwph 



dwph 



where 



dPeli'Pf.9f) 

dwph 

rffmag.l(y/,6'/) 

dwph 

rf-Pmag,2(tp/,6'/) 
dWph 

d -Pinterforcnco(<y2/ , 9f) 



72 2 


WphE^ 


27rc5 


vn? ki 


ry2 2 


Wpb Ei 


27rc5 


m?ki 


72 „2 


WphE^ 


27rc5 


rn? ki 


72 2 


■WphE, 



Pe.(fc.fc/) '^-f''^'"^^ h.C. 

d cos 9f 
dp*iih,kf,nf) 



Pmag,! \ki, kf 



Pmag,2(fcj,fc/)- 



d cos9f 

dp*ynagAh,kf,nf) 



h.e. 



d cost 



7 



dw 



ph 



27r c^ m"^ ki 

Pmag,l i"^i, f^f) 
Pmag,2 ih,kf) 



Pc\{ki,kf 



d (Pmag,l(^*' kf,^f) +Pmag,2(fc*> %> ^/)) 



d cos 9f 

d {pl\{ki, fc/, »/) +p;;,ag,2(fc»: kf.nf)) 

d cos 0/ 

cf ipl\{k.„ fc/, %) +p;,ag i(fc,, fc/, %)) 



d COS0 



h.e. 



/ 



(51) 



For clarity of further analysis, we call dPc\ as electric component of emission (or electric emission), dPmag,! as 
magnetic component of emission (or magnetic emission), dPmag,2 as correction of magnetic component of emission 
(or correction of magnetic emission), d Pintorforoncc as interference component of emission. Sometimes, we shall omit 
variables (/?/, 9f in brackets of these functions. 

For description of the bremsstrahlung which accompanies collisions of protons off nuclei, we shall consider in this 
paper only normalized cross-section as 



d'^a 



dwph d cos 9f 



= NoWpY, ■ <p{ki,kf) 



^dp *{ki,kf,Vtf) 
d cos 9f 



h.e. 



(52) 



where A'o is normalization factor (determined by normalization of the calculated curve of the full bremsstrahlung 
spectrum on 1 point of experimental data), and in calculations of matrix elements we use boundary condition of 
elastic scattering for the wave function of the proton-nucleus system in the state before emission of photon. 



III. RESULTS 



Let us estimate the bremsstrahlung probability accompanying the proton-decay. I calculate the bremsstrahlung 
probability by eq. ()49p . The potential of interaction between the proton and the daughter nucleus is defined in 
eqs. (26)-(27) with parameters calculated by eqs. (28)-(29) in [73- The wave functions of the decaying system in 
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the states before and after the photon emission are calculated concerning such potential in the spherically symmetric 
approximation. The boundary conditions and normalization are used in form of (B.1)-(B.9) in |75j . To choose the 
convenient proton-emitters for calculations and analysis, one can use systematics presented in Ref. [J| (see Table II 
in the cited paper). In particular, in [75| the "'^^^Ta, ^^"'^Re, ^^''Ir and ^®^Bi nuclei decaying from the 2si/2 state (at 
k = 0), the 53^156, 5pCs57 nuclei decaying from the 1^5/2 state and the g^^Tmry, yf^Luso nuclei decaying from the 
0/iii/2 state (at k ^ 0) were selected. In this paper I shall analyze only one nucleus g|^Tm77 at li ^ (as calculations 
for this nucleus are essentially more difficult than for nuclei at h — 0), with a main emphasis to study new physical 
effects in frameworks of the proposed model (assuming that such studied effects should be similar for other nuclei). 
For the IfTuijr nucleus we have h = 5, If = 4, Q = 1.140 MeV [zl. 



A. Electrical, magnetic emissions and angular distributions 

At first, let us clarify how much the magnetic emission is visible on the background of the full bremsstrahlung 
spectrum (to understand if there is a sense to study it, at all). The result of calculations of the bremsstrahlung 
probabilities during proton decay of -'^^^Tm (at the chosen angle 6 = 90° between the directions of the proton motion 
(with its possible tunneling) and the photon emission) are presented in Fig. [TJ The electric and magnetic components 
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FIG. 1: (Color online) The full bremsstrahlung spectrum, electric and magnetic components of emission defined by eq. (|5ip 
(at 6 = 90°): (a) the full spectrum (full blue line), electric component dPd (red dashed line) and magnetic component dPmag.i 
(green dash-dotted line), (b) ratio of the components to the full spectrum (full blue line is for dPd/dPfuii, red dashed line for 
dPina.g,i/dPtuu)- One can see that the magnetic emission gives contribution about 28 percents inside energy region 50-300 keV. 



are included also on these figures. One can see that the magnetic emission is smaller than electric one. But it gives 
contribution about 28 percents into the full spectrum (see Fig. [T] (b)), i.e. it is not so small to be neglected and it 
should be taken into account in further calculations of the bremsstrahlung spectra during nuclear decays with emission 
of charged fragments with non-zero spin. However, the magnetic component suppresses the full emission probability: 
according to Fig. [T] (b) (see blue solid line), inclusion of the magnetic component into calculations is determined by 
Pei/Pfuii — 1.14, which is larger unity. This effect of suppressing of the total emission can be explained by a presence of 
not small destructive interference between the electric and magnetic components inside whole studied energy region. 
According to Fig. [T] (b) , ratios of the electric and magnetic components to full spectrum are not changed in dependence 
on the energy of the emitted photon. As we find, the correction of the magnetic component dPmag,2 is smaller than 
the electric and magnetic components by 10^ times (so we shall neglect by such a contribution in further analysis). 

Now we shall analyze how the magnetic emission is changed on the 9 angle between the outgoing proton and emitted 
photon. In particular, let us find if there are some values of such angle, where the magnetic emission increases strongly 
relatively electric one. In Fig. [2] the angular distributions of the electric and magnetic emissions during the proton 
decay of ^"^^Tm are shown. One can see that the electric and magnetic components increase proportionally (similarly) 
with increasing of the 9 angle. From Tabl.|T]it follows that there is no any angular value, where the magnetic emission 
increases essentially relatively electric one. 
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FIG. 2: (Color online) The angular distributions of the bremsstrahlung emission during proton decay of the ^''^Tm nucleus: 
(a) the electric component of emission, dPd, calculated at different energies of the emitted photons; (b) the electric component 
dPei (full blue line) and magnetic component dPmag.i (red dashed line) for the chosen photon energy 200 keV. One can see that 
both spectra increase proportionally (similarly) with increasing of the angle. 



Angle 

e 


Emission probability | 


Electric 
component dPd 


Magnetic 
component dPmag,i 


dPn.ag,l/d^cl 


10° 
20° 
30° 
40° 
50° 
60° 
70° 
80° 
90° 


1.704 X 10"'* 
2.580 X 10"" 
1.192 X 10"" 
3.329 X 10"" 
6.952 X 10"" 
1.188 X 10"" 
1.727 X 10"" 
2.158 X 10"" 
2.319 X 10"" 


4.198 X 10"'° 
6.357 X 10"" 
2.940 X 10"" 
8.212 X 10"" 
1.716 X 10"" 
2.939 X 10"" 
4.281 X 10"" 
5.361 X 10"" 
5.779 X 10"" 


0.24630 
0.24636 
0.24647 
0.24665 
0.24692 
0.24730 
0.24779 
0.24841 
0.24916 



TABLE I: Electric and magnetic components of emission in dependence on the 9 angle between directions of outgoing proton 
and emitted photon at 200 keV of photon energy. One can see that ratio of magnetic emission on electric one is not changed 
practically inside whole angular region. 

B. How are the electric and magnetic emissions changed in dependence on distance between the proton and 

the daughter nucleus? 



Usually, authors of papers on the bremsstrahlung, which accompanies different types of collisions of particle between 
themselves and with nuclei, decays and fission of nuclei, calculate the spectra on the basis of integration over all 
space coordinates. In the relativistic models of collisions of nucleons off nucleons and nuclei (at intermediate energies) 
calculations are preformed in impulse representation mainly. Such approaches miss information on how much intensive 
emission is in dependence on distance between centers-of-masses of two studied objects. However, it is natural to 
think that photons are emitted with different intensity in dependence on such a distance. One can suppose that 
electric and magnetic photons are emitted by different ways. We put such questions: 

1. Can the magnetic emission be stronger than electric one in some space region? 

2. How are the electric and magnetic emissions changed in dependence on distance between the proton and nucleus? 

3. How much strong are the electric and magnetic emissions from the tunneling region? Is there principal difference 
between these types of emission from the tunneling region in comparison with the external one? 

In order to perform such an investigation, we shall define the probability of emission of the bremsstrahlung photons 
from the selected space region. In the presented formalism the dependence of emission on the distance is determined 
by the radial integrals Ji(^i, lf,n), J2{h, If, n) and Jaih, If, n) in eq. (j39|) and (j4ip . where integration if performed over 
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full space region. So, to obtain emission from arbitrary selected interval r G [ri,r2], we shall consider the following 
integral: 



Jm{lijf,n;ri,r2) 



fmir) dr, 



(53) 



where m — 1,2,3 and fm{f) is integrant function of the corresponding radial integral Jm{h, '/j '^)i defined in eq. p9p 
or (j4ip . In particular, Jm{hilf,n]ri,r2) transform to Jm{li,lfTn) at ri — > and r2 — )• +oo. Now, for the emission 
from enough small interval Ar near studied distance r we obtain: 



r+Ar 



Jmikjf,n;r,r + Ar) = 



frn{r') dr'. 



(54) 



From here we define amplitude of emission in dependence on the distance r on the basis of such a radial function: 



JjnikJf,ri;r,r + Ar) 



r+Ar 



r+Ar 



Ar 



Jm{k,lf,n\r) = lim '■ — '—^ = lim -t— 

Ar— >0 

= Ui{r) lim -^ Ar = f„^{r) 



r-i-O Ar J Ar^O Ar J 



(55) 



1 



Ar^O At 



After this, matrix elements and probability of emission with included dependence on the distance r can be defined as 
before, where we shall use Jmih, If, n; r) instead of the radial integrals Jrn{h, If in). For denotation of new character- 
istics with dependence on the distance r we shall include variable r inside brackets. 

The magnetic component dP^ing i(r) on the background of the electric one dPci{r) in dependence on the distance 
r is shown in Fig. [3l One can see that behaviors of both functions are similar: they oscillate in the external region 
(having maxima and minima at similar space locations), while they have monotonous shapes with one possible well 
in the tunneling region. In general, the magnetic emission suppresses the full emission inside whole space region. The 
emission from the internal region up to the barrier is the smallest, and from the external region — the strongest. 
Behavior of the correction of the magnetic component (iPmag,2('') on the background of the electric one dPei{r) in 
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FIG. 3: (Color online) The magnetic component dPniag,i(r) and the electric component dPci{r) in dependence on distance 
r between centers-of-masses of the proton and daughter nucleus at 200 keV of the emitted photon energy (at 6 = 90°): (a) 
the magnetic component dPmag,i(r) (red dashed line) and the electric component dPci(r) (full blue line) inside the space 
region up to 250 fm. One can see that both functions oscillate similarly in the external region outside the barrier, while they 
are essentially smaller inside the tunneling region; (b) the magnetic component dPinag,i(r) (red dashed line) and the electric 
component dPci{r) (full blue line) inside the tunneling region (up to 80 fm). One can see that both functions have monotonous 
behavior (with possible one well, and without any oscillation) in this region. After crossing from the barrier region into the 
external one the first oscillation is appeared with further peak sharply increased (this demonstrates more intensive emission 
from the external region in comparison on the tunneling region). One can see also that after crossing from the barrier region 
into the internal one (near 12 fm) strong decreasing of both functions is appeared (with oscillations) — this points on the 
extremely smaller bremsstrahlung emission from the space region of nucleus; (c) ratio of the magnetic component to the electric 
one, dPmag,i(r)/dPci(r) (full green line). One can see that this characteristic is not changed inside whole studied region of 
distances, it is the same in the tunneling and external regions. 
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dependence on distance r is shown in Fig. |4l In general, this function is essentially smaller. In the tunneling region 
it increases monotonously, in contrast to the electric and magnetic components (see Fig. |3](c)). This causes a sharp 
peak of the function dPma.g,2{r) / dPc\{r) close to the external boundary of the barrier (external turning point) shown 
in Fig. 131(b). This peak could be interesting for further research, as it corresponds to the external space boundary of 
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FIG. 4: (Color online) Correction of the magnetic component of emission, dPmag,2(r), in dependence on the distance r between 
center-of-masses of the outgoing proton and the daughter nucleus at 200 keV of energy of the emitted photon (at 9 = 90°): (a) 
the correction of the magnetic component dPmag,2(»') x 10^ (red dashed line) and the electric component dPci{r) (full blue line) 
in region up to 250 fm. One can see that in the external region outside the barrier both functions oscillate similarly, while in 
the tunneling region (up to 80 fm) they are essentially smaller; (b) ratio of the correction of the magnetic component to the 
electric one, dPina.s,2{r) / dPci{r) . One can see that there is a sharp peak close to 80 fm (that corresponds to the external turning 
point); (c) the correction of the magnetic component dPina.g,2{r) x 10^ (red dashed line) and the electric component dPei{r) 
(full blue line) in the tunneling region. One can see that in this region these two functions has principally different behavior. 
By this difference one can explain the presence of the peak in the previous figure (b). 

the barrier. But, unfortunately, this peak is extremely small (in comparison with the full spectrum) for any reasonable 
searches of its experimental measurements. 



C. Spectra of the emitted soft photons 



From point of view of theory, it can be interesting to know what happens with the bremsstrahlung spectrum at 
limit of energy of the emitted photons to zero. In particular, let us analyze if this spectrum increases infinitely or 
tends to definite finite value, and which limit is in that case. 

For low energies of photons (i.e. for soft photons) two prevailing approaches are known: the first approach is started 
from the early work [8l| of Low and it is based on application of soft-photon theorem to all nuclear bremsstrahlung 
processes, the second one is based on application of the approximation of Feshbach and Yennie :82], which is more 
effective near resonances (see [Ij] for analysis). However, as it was noted in [ij] (see p. 376), there is another way 
of development of the bremsstrahlung theory, i.e. potential one, to which our model can be referred. According to 
theory of QED, divergence in calculation of the matrix element is appeared at limit of photon energy to zero (infrared 
catastrophe, see p. 258-273 in [60]; P- 194-200 in ^3|; P- 194, 225, 231 in [SOl)- However, we obtain the convergent 
integrals and the finite probability of the bremsstrahlung emission in our approach. In particular, let us consider the 
first integral in eqs. (|39|) for n = at limit uiph — >■ 0: 



lim Ji{li,lf,n ^ 0) = lim 



-Ro=i/fcph 

dRi{r,li) 



+ 00 



dr 



R*fih^'^)Jo{kphr) r^dr 



lim 

Uph^O 



dRi (r, k 



dr 



-Ro=i/fcph 



(56) 

At Wph — >■ we have jo(k-phr) — sin(fcphr)/(fcphr) — >■ 1 (fcph = Wph/c). So, one can see that the first item converges 
(according to the chosen boundary conditions, X/('') = at r = 0, where Rf{r) = X/('^)/'' IZSil)- The second item 
does not include small energies of photon (fcph > l/i?o) and, therefore, it is standard integral in our calculations of 
the spectra of not-soft photons, i.e. it converges also. The same result can be obtained at arbitrary chosen n and 
for J2{li,lf,n), J {li,lf,n). On such a basis, according to eqs. (j36l) . (|38)) and (|40)) . all matrix elements Pei, Pmag,i 
and Pmag,2 (and the angular matrix elements) converge at arbitrary values of quantum numbers li, If. According to 
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eq. (|49t . we obtain: 






lim 
Uph^o dwph 



h.e. 



0. 



(57) 



d cos9f 
Our calculations at energy of the emitted photons up to 2.5 keV are shown in Fig. [5l One can see that at decreasing 



>; 10"^- 




146^. 

Tm 


■a 
> 

^ 10"'- 
■a 

D. 

■a 

10"°- 


f 


(a) 



0.5 1.0 1.5 2.0 

Photon energy, E (keV) 



i 10^- 




146^ 

Tm 




^f^:^^ -^'^^:^~~ 


10"'- 


t - - ; 


^:'^ 


- 10"'- 


I'f -.-. 






(b) 


10"'- 


/ ,.-"""■" "^"""""•-.- 


-. 



0.5 1.0 1.5 2.0 

Photon energy, E CfceV9 



FIG. 5: (Color online) The bremsstrahlung spectra for near-zero energy of the emitted photons (up to 2.5 keV): (a) full spectrum 
(full blue line), electric component dP^i (red dashed line) and magnetic component dPmag.i (green dash-dotted line) at 6 = 90°; 
(b) full spectrum in dependence on the Q angle (full blue line for Q = 90° , orange dash-dotted line for Q — 75° , red dashed line 
for 9 — 60°, olive short dotted line for 6 — 45°, wine short dash-dotted line for 6 = 30°, violent short dashed line for 9 = 15°). 

of the photon energy the bremsstrahlung probability increases slowly up to finite maximum, and then it decreases 
monotonously. According to our estimations, the probability has finite maximum at energy of the emitted photon 
smaller 1.5 keV. So, there is no the infrared catastrophe in our approach"^. 

D. Spectra in collisions of protons off nuclei at intermediate energies 

In finishing, I shall shortly demonstrate applicability of the proposed model and calculations for description of 
experimental spectra of the bremsstrahlung during collisions of protons off nuclei at intermediate incident energies of 
protons. I calculate the normalized cross-sections by eq. ([5^. use the same form of the proton- nucleus potential and 
parameters (defined as for the problem of proton-decay studied above)^. 

In Fig. [5] (a) one can see that our approach can describe enough well experimental data for p +^ Be in energy region 
from 20 MeV to 120 MeV in comparison with results obtained by Nakayama and Bertsch in [SJ] and calculations 
performed by Nakayama in [6l|. In next figure (b) we compare our calculations for p +^'^ C with experimental 
data [85j] and results obtained by Remington, Blann and Bertsch in [86|. Such comparison shows: in energy region of 
the emitted photons up to 90 MeV our full spectrum (see solid blue line) is enough close to experimental data and 
calculations obtained using master equation and semiclassical bremsstrahlung formula (see wine dash double-dotted 
line) , the semiclassical cross-sections multiplied by factor 2 for meson exchange (see orange short dotted line) in [8al ■ 
But for hard photons with energy from 90 to 120 MeV we achieve better agreement with experimental data than 
results of [8a] . Comparison of our results with quantum calculations performed in [8^ (see navy short dashed line 
in that figure) indicates on absolute applicability (and availability) of the quantum approach in description of the 
emitted photons of high energy in collisions of protons off nuclei. At the same time, such an approach allows to deeper 
study quantum properties (such as non-locality, for example) of the considered colliding process. In last figure (c) 



It is interesting to note that such a proposed definition of probabihty, eq. I I49I I , allows to describe enough well experimental data of the 
bremsstrahlung emission during a-decay without any normalization of the calculated spectra on experiment (see Fig. 1 in [331). 
A key problem in obtaining the reliable bremsstrahlung spectra is difliculty to achieve stability in calculations of the matrix elements. 
Also I suppose that this is the main reason why a main idea of the proposed potential approach was not developed essentially for 
calculations of the bremsstrahlung spectra at intermediate energies which accompany different kinds of nuclear processes. In order to 
achieve the stability, I applied the approach presented in Appendix in |49l | inside the radial region from i?as to fJmax- For simplicity of 
analysis, I used the same values for these two parameters: Has = 0. 9 X {Rn + 7cir), iJmax is chosen so large when the spectrum is not 
changed after its variations, Rji and ajj are potential parameters defined in eqs. (29) of | 75| . 
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FIG. 6: (Color online) The proton nucleus bremsstrahlung probability rates in the laboratory system at the incident energy 
31ab ~ 140 MeV (in our calculations we use photon emission angle 6 — 90°): (a) Comparison for p+^Be between the calculations 
by our model (blue solid line is for full spectrum, green dash-dotted line for electric contribution, red dashed line for magnetic 
contribution), calculations from (Nakayama 1986: [SJ], wine short-dashed line), calculations from (Nakayama 1989: 61], navy 
dash double-dotted line) and experimental data (Edington 1966: [89]); (b, c) Comparison for p+^^ C and p-l-^"* Pb between the 
calculations by our model (blue solid line is for full spectrum, green dash-dotted line for electric contribution, red dashed line 
for magnetic contribution), calculations by Remington, Blann and Bertsch in (Remington 1987: [86|], wine dash double-dotted 
line is for calculations by master equation using the semiclassical bremsstrahlung cross sections, orange short dotted line for 
semiclassical cross sections multiplied by 2 for meson exchange, and navy short dashed line for quantum bremsstrahlung cross 
sections), and experimental data (Edington 1966: [Sa]'). 



similar comparison is performed for p -|--"^ Pb. On all figure we add our calculations for the magnetic and electric 
bremsstrahlung emission. 

In Fig. |7| we present our calculations of the bremsstrahlung cross-sections for collisions p +^ C, p +^'^ Cu and 
p_l_i07^g jj-^ comparison with experimental data ,87,] at the incident proton energy Tiab = 72 MeV. Here, we show the 
full spectrum calculated by eq. (|5^ and the corrected spectrum obtained by eq. (|5^ with division on fc/ (according 
to formula (13) of cross-section defined in 56]). Comparison with quantum calculations performed by Kopitin, 
Dolgopolov, Churakova and Kornev in ^56] (see Fig. 1 in the cited paper) shows more stable calculations in our 
approach. In addition, this answers on assumption putted in [56| that quantum approach (with included nuclear 
component of potential) is absolutely able to describe well experimental data of the bremsstrahlung during proton- 
nucleus collisions. 
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FIG. 7: (Color online) The proton nucleus bremsstrahlung probability rate in the laboratory system at the incident energy 
Tlab = 72 MeV and photon emission angle 6 — 90°: Comparison for p +^^ C (a), p -t-®** Cu (b) and p -f i"'' Ag (c) between the 
full cross-section calculated by eq. (|52p (wine dash double-dotted line) , the corrected cross-section obtained by eq. (|52p with 
division on kf (blue solid line) and experimental data (Kwato 1988: [83] )• We add the electric component (green dash-doted 
line) and magnetic component (red dashed line) to all figures. 



Comparing results of calculations obtained for p+-^'^ C, p+^'^ Cu and p-|-^°^ Ag at the incident energy Tiab = 72 MeV, 
one can find worse agreement between theory and experiment for the ^^Cu nucleus. This situation looks to be enough 
strange, as all measurements were made by the same group of experimentalists (and it is difficult to expect that the 
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cross section for the ^^Cu nucleus is the experimental error). For evidence, let us consider all these experimental data 
in one figure (see Fig. |S]). One can see that data for ^''Cu are located lower than the data for ^^C and ^°'^Ag. At the 
same time, the data for ^''Cu are decreased more slowly with increasing of the photon energy than the data for ^^C 
and ^°^Ag. In particular, one can expect that further continuation of all these data for higher photon energies have 
to lead to their intersection at one point (or these are evident deviations from monotonous decreasing tends of the 
spectra), that has never been observed before. 
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FIG. 8: (Color online) The experimental data (Kwato 1988: ^) for p +^^ C, p -h"" Cu and p +'^°'' Ag at Tiab = 72 MeV: (a) 
One can see that data for ^''Cu are located lower than the data for ^^C and ^"^Ag. At the same time, the data for "^Cu are 
decreased more slowly with increasing of the photon energy than the data for ^^C and ^°^Ag. (b) The comparison between 
experimental data reinforced by calculations of the full cross-sections (blue solid lines in Fig. [71 = 90°): inclusion of the 
calculated curves, describing general tendency of the spectra, only reinforces difference in behavior between experimental data, 
indicated on the previous figure (a). 

Such a picture disagrees with early observed general tendency of the bremsstrahlung spectra in nuclear processes. 
This logics explains why the spectrum is decreased more strongly (with increasing of the photon energies), if this 
spectrum is lower. Such a tendency is based on correspondence between shape of the barrier with tunneling length of 
the emitted fragment: the energy of proton is lower, the length of tunneling is larger, the total emission of photons 
is less intensive (because of it is less intensive for photons emitted from tunneling region, than from above-barrier 
regions). We demonstrated this tendency on the example of two spectra in the a-decay of the ^^^Po and ^^^Ra nuclei 
(see Fig. 3 and explanations in [33]). 

But, if this difference between experimental data is supported by future measurements, then such a result would 
be very interesting. This will be a direct indication on influence on the spectra of some other hidden characteristics 
of the proton-nucleus system, which are not included into current calculations. This will the indicate on a presence of 
new aspects in the bremsstrahlung spectra. One can assume that structure of the proton-nucleus system, dynamics 
of its nucleons, other early not studied properties can be important in such new developments^. 

E. Role of the multipolar components in the angular analysis 

The first calculations of the multipolar components of the bremsstrahlung emission of higher order in tasks of 
nuclear decays were obtained by Tkalya in |25l . l26j . Studying emission of the bremsstrahlung photons during a-decay 
of the ^^^Ra, ^^°Po and ^^''Po nuclei, he shown that the multipolar term E2 is essentially smaller in comparison with 
El (see Fig. 1 in [2g, ratio between contributions p^^/p^^ is about 50-1000 for the photon energy range up to 
900 keV). There are also estimations obtained by Kurgalin, Chuvilsky and Churakova for the multipolar term E2 
of the emitted photons in a-decay of ^-'^°Po [Sg: according to their calculations, contribution of the E2 multipole is 
smaller than El by 50-500 times for the photon energies up to 800 keV. We studied this question also and found 
the multipolar terms E2 and M2 to be very small. Authors of paper [35| investigated the dipole and quadrupole 
contributions in the semiclassical consideration to the bremsstrahlung probability in a-decay, studied interference 



^ For example, in the problem of the bremsstrahlung emission accompanying ternary fission of ^^^Cf (where this nucleus is separated 
on the a-particle and two heavy fragments) we shown that dynamics of relative motion of all participated fragments, and geometry of 
nuclear separation have strong infiuence on the bremsstrahlung spectrum [52| . 
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between such contributions^. I have not found any information about other attempts to estimate the E2 multipolar 
term and the multipoles of higher order, which could be obtained up to now. By such reasons, calculations of the 
bremsstrahlung spectra in the multipolar approach usually are performed on the basis of the first multipolar term, 
which gives the prevailing contribution into the full spectrum (usually 4-5 first digits of the calculated probability are 
stable in our approach, as minimum). 

Also it is more difficult to obtain reliable estimations of the multipolar terms of higher order because of essentially 
smaller convergence of their calculations. This is real practical difficulty (which can be alienated from many researchers 
trying to obtain the multipolar terms of higher order). Indications on difficulty of such problems and perspective of 
their solution I find in papers of authors, who calculated the bremsstrahlung spectra in different nuclear tasks with 
realistic potentials (for example, see |13l. Ilil. ISq ) . 

In order to understand more clearly, how the angular bremsstrahlung probability is changed in dependence on 
quantum numbers h, If and Zph (which defines the multipolar term), we rewrite formulas separating components 
which describe this angular dependence. This information is completely included in the differential matrix elements: 



dpf' 



'ph/j 






dp? 



smOde 



g pi"/ 1 Y^ 



P, 



\fj.-fj.' 



M' = ±1 



\S, jn r^' pl™*"^'l + r^' pl™'-/^' 






u' p|m,-A''l 
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P, 



Im-m'I 
lph-1 



(58) 
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'r'{^8<.-i-^^^.:.+i}> 



(59) 
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(61) 



Expansion in 35| and the multipolar expansion in the given paper have different basis and sense. In 13511 dipole and quadrupole 
contributions are defined as the first term (at Ij = 1) and the second term (at If = 2) of expansion of wave function (/3 j (r) of the a- 
nucleus system in the state after emission of photon (see eqs. (B1)-(B4) in [35l| ). at representation of the effective charge for two-charged 
nuclear system (see eqs. (A1)-(A4) in [331). The multipolar approach in this paper is based on the standard multipolar expansion of 
the wave function of photon II28I I. 
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ph 



I ^iii/iph,iph-i ■ "^v *' /' ph ^)' 



(62) 



^ph 
tph 



2/n7TT ^'"'/'ph.'ph+i ' "^ ^^" ^■''' ^p'' "^ "^''' 



Here, c^ ... Cg and cj ... cg are not dependent on the 9 angle. Function Si-^q is defined as Si-^q = at /^ = and 
Sii^o = 1 at Ij ^ 0. Formulas for the first some values of li and If are presented in Appendix [C] On the basis of these 
formulas we conclude the following. 



1. Numbers li and If determine basic shape of the angular distribution of the bremsstrahlung probability, number 
Zph determines oscillations in this shape: 

(a) Number of oscillations of this shape is minimal at Zph = 1 and increases at increasing of /ph. 

(b) Ci ■ ■ ■ Cg and C7 . . . cg are weights of oscillations at each chosen Zph. As integrals Ji, J2 are decreased at 
increasing of Zph (at fixed Wph), so each matrix element with next value of /ph gives own new contribution 
into base shape of the probability distribution with smaller intensity, but larger number of oscillations. 



,-m'I 



I™/ 1 



at 



2. If polynomials P;™\ ^ at some chosen li or polynomials P; at chosen // in eqs. ((58|) (polynomials P; 

some chosen li or polynomials P^''"^' at the chosen If in eqs. ([5^ ) equal to zero for some values of the 9 angle, 
than the differential matrix elements in eqs. ((58)) (in eqs. ([59)) ) equal to zero at any value of /ph for this 9 angle. 

The angular contributions of the electric component dPd of the bremsstrahlung emission during the proton decay of 
the ^^^Tm nucleus for the first three multipoles are presented in Fig. [9] In figure (a) one can see that the second and 
third multipolar contributions (at /ph = 2 and /ph = 3, = 90°) are smaller on 5-7 orders of magnitude in comparison 
with the first one (at Ipi^ — 1, 9 — 90°). The angular distributions of these multipolar contributions are shown in next 
figures (b, c) for /ph = 2 and /ph = 3. In particular, one can see that for smaller values of the 9 angle the emission is 
more intensive at increasing of the multipolar order /ph (at the same fixed h and // for ^ Tm). 
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FIG. 9: (Color online) Contributions of the electric component dPd of the bremsstrahlung emission for proton decay of the 
^*®Tm nucleus for the first three multipoles {Iph = 1,2,3). (a) The spectra at 6 = 90°: one can see that the first contribution 
at Iph = 1 (blue solid line) is essentially larger in comparison with contributions at /ph = 2 (green dashed line) and /ph = 3 (red 
dash-dotted line), i. e. the first multipolar contribution is prevailing inside whole energy region of the emitted photons, (b) 
The multipolar contribution at /ph = 2 in dependence on the 6 angle: one additional extremum can appear in each curve inside 
the angular region from up to 90°, but it is practically smoothed (at current computer accuracy of calculations). However, 
at small values of 9 each curve is increased more sharply in comparison with the angular spectra at Iph = 1 (see Fig. [2] (a)). 
(c) The multipolar contribution at /ph = 3 in dependence on the 9 angle: appearance of else one new extremum in each curve 
forms one new oscillation. There is displacement of maximum and minimum of each spectrum in direction of larger values of 
9 with increasing of energy of the emitted photons. 
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IV. CONCLUSIONS 

The new model of the bremsstrahlung emission which accompanies proton decay and coUisions of protons off nuclei 
in the energy region from the lowest up to intermediate, has been developed. This model includes spin formalism, 
potential approach for description of interaction between protons and nuclei, and operator of emission includes the 
component of magnetic emission (defined on the basis of Pauli equation). In the problem of the bremsstrahlung 
during the proton decay in the first time a role of the magnetic emission is studied using such a model. For such 
investigations the ^^^Tm nucleus is chosen. We obtain the following: 

1. Inside energy region from 50 up to 300 keV the magnetic emission gives contribution about 28 percents into 
the full spectrum (see Fig. [T|), i.e. it is not so small and should be taken into account in further estimations of 
spectra of the bremsstrahlung emission during nuclear decays with emission of the charged fragments with non- 
zero spin. However, the magnetic component suppresses the full emission probability: inclusion of the magnetic 
component into calculations is determined by Pd/Pfuii — 1-14, which is larger unity. This effect of suppressing 
of the total emission can be explained by a presence of not small destructive interference between the electric 
and magnetic components inside whole studied energy region. Ratios of the electric and magnetic components 
to full spectrum are not changed in dependence on the energy of the emitted photon. The correction of the 
magnetic component dPmag,2 is smaller than the electric and magnetic components by 10^ times. 

2. With increasing of the 6 angle between directions of the outgoing proton and emitted photon the electric and 
magnetic components increase proportionally (see Fig. [5]), but ratio between them is not changed (see TabLH]). 
So, there is no some angular value, where the magnetic emission increases essentially relatively electric one. 

3. The magnetic component (i-Pmag,i('') is dependent on distance r between centers-of-masses of the proton and 
daughter nucleus similarly as the electric component <iPmag,i(?') (ratio between such two components is not 
changed inside region from 5 fm up to 250 fm). In the external region both components oscillate (having 
maxima and minima at similar space locations), while in the tunneling region they have monotonous shapes 
with one possible well (see Fig. [3]). In general, the magnetic emission suppresses the full emission inside whole 
space region. The emission from the internal region up to the barrier is the smallest, and from the external 
region — the strongest. 

4. The correction of the magnetic component dPmag,2(?') is essentially smaller than the electric one dPci{r) in 
dependence on distance r (see Fig. 3]). In the tunneling region it increases monotonously, in contrast to the 
electric and magnetic components. This causes a sharp peak of the function dPma.g^2{i') / dPc\{r) close to the 
external boundary of the barrier (near 80 fm). 

5. At decreasing of the photon energy up to zero, the bremsstrahlung probability increases slowly up to finite 
maximum (at energy of the emitted photon less 1.5 keV), and then it monotonously decreases to zero (see 
Fig. [5]). The angular distribution of the probabilities of the bremsstrahlung emission at such small energies 
looks like the angular distributions inside the energy region from 50 to 350 keV studied above. We show that 
there is no the infrared catastrophe in our approach. 

It is demonstrated that the model is able to describe enough well experimental data of the bremsstrahlung emission 
which accompanies collisions of protons off the ^C, ^^Cu and ^''^Ag nuclei at the incident energy Tiab — 72 MeV (at 
the photon energy up to 60 MeV), the ^Be, ^^C and ^''^Pb nuclei at the incident energy Tiab — 140 MeV (at the 
photon energy up to 120 MeV). 

Appendix A: Linear and circular polarizations of the photon emitted 

Rewrite vectors of linear polarization e'"' through vectors of circular polarization ^^ with opposite directions of 
rotation (see Ref. (zll, (2.39), p. 42): 

e_i = i=(eW-ze(2)), e+i = --^(eW+.e(2)), ^, ^ e^^) , (Al) 

where 

h±=T^^, h^i+h+i^-iy/2, E e(")'* = /i-id + /i+iCi- (A2) 



V2 



a=l,2 
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We have (in Coulomb gauge at e'^^^ = 0) 



1 
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We shall find also multiplications of vectors ^j.^. From eq. (jATj) we obtain: 
From here we find: 
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Appendix B: Angular integrals Ie, Im and / 

We shall calculate the integrals in eqs. (|39p and (|4T|) : 

/m (/^, ;/, ^ph, h,^^) =1 >i;„,, (n,.) T,, i,^m. (n,) T,;,^ ,^,^^ ^(n, 
lEih,lf,lphJi,l2,tJ-) = J Yi*.mfinr)Ti^i,^rn,{n,)T'l^^i^^^{nr) dn, 
I{h,lfjph,n,ti) = 1^ / Y,*,„^(nr)r(,™,(nr)T;*^,^„^^(nr)drj. 

Substituting the function Tji^m{n,.) defined by eq. ([27]), we obtain (at £_q — 0): 

Im {h,lf,lph,li,f^) = ^ {li,l,li\mi- fj.',fi',mi) {lph,l,lph\fJ-- ^J■',^J■',^J') >^ 

X JYi].„,{n,) ■ yi„,„,_^,(nr) • r4,^_^,(n,) dn, 

lE{k,lf,lph,h,l2,fJ') = ^ (?i,l,?., I mi - ^',^',mi) (?2,l,^ph I M-/^':M',m) X 

m'=±i 



(Bl) 



(B2) 



I{k,lf,lpi,,n,n) = (n,l,/ph|0,^,Ai) X / y;*„^(nr)11,™,(nr)r„*o(nr) dr2. (B3) 

Here, we have taken orthogonality of vectors ^+i into account. In these formulas we shall find angular integral: 
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where Pj"'(cos0) are associated Legandre's polynomials, and we obtain conditions: 

for integrals Im,Ie ■ fJ. = mi — mf, n>\fi - fi'\ = |m,: — m/ + /i'|, ^ = ±1, 
for integral / : rrii = irif. 

Using formula (|B4p . we calculate integrals (jB2p and (jB3P : 

7r 

TT 



where 



^mim/A' _ /•_-[^y/+i2+mi-A'' jlf+h+l2 + \mf\ + \m,-tj.'\ + \m^-mf~ti,'\ ^ 
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We define differential functions on the integrals (jB6[) with angular dependence as 
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m'=±i 



'^u™-™, y c,77^;, •/r;r^(e), (BIO) 

SmOdO ' l■^l■fl■phn l^Lfn 

Appendix C: Differential matrix elements for the fist U and If 

We write calculations for the fist some values of U and //, at arbitrary /ph: 
1. k ^0,lf^ 0: 



2. /, = 0, // - 1 



rfP/*^ 



Af 



sin6'd6' 



s^ - ~ sin^ ^ 4'' . P,'^7'''', m. = 0, m/ = ±1, 



— sm I 



E {^5' <-"i'' - <' ^ir/i'' }' ™^ = 0, m, = ±1, 



— ^ = cj ■ COS0 P,7 , m,; = 0, m/ = 0, 

sm 9 do p'' 

c^pf „ r -I 
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3. h^OAf 



2: 






sm0d0 

~M 
'phM 



dvf 



sin Ode 
dpf 

smOde 



- 3 sin^ 61 cos 61 ^ c^' • P, 
3 8102 6*008 61 ^ |c^' 

At'=±l 

1 

2 



ph 



^ph — 1 6 /ph + 



771, = 0, mf ~ ±1, 



5''}. 



0, TO/ = ±1, 



(C3) 



TOi ~ 0, TO,/ = 0, 

\ (3 cos^ - 1) {c8 P,"^^_i - eg pO,^+i}, m, = 0, TO/ = 0. 



4. L^l /f = 1: 



dp^ 



^phM 



8in 61 dd 

^phM 



C08 6'- ^ i^-Sm.fiic'^ SVCL^ + 5ra,,±l Sm,fi' c"^" ~ C2 P] 



M' plTli-^'l 



•P 



Im-p' 



'^Pf 



8in d0 



= CO8 • 2, { ^niifi 3 8in 



^2 /) V 



3mi,±l 



C08 6* • ^ <^ (5„,_o 3 sin^ c'^ + S,n,,±i 






-'nii^ '^3 



^rni/j.' '-'4 I *-'6 2 



P 



Im-p' 



P 



iph-i 



iph+1 



-'^'phA' 

8in d0 
8in d0 



C,-{p\r\fPl^, TO.=TO/=0,±1, 

(Pl"'l)'{c8 P0^_, - eg P0,^ + ,}, TO. = TO/ = 0, ±1. 



-li - TO/ I 



= 1, 



(C4) 

(C5) 
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